Square vortex lattice at anomalously low magnetic fields 
in electron-doped Nd x 85 Ce .i5CuO4 
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We report here on the first direct observations of the vortex lattice in the bulk of electron-doped 
Ndi.85Ceo.i5Cu04 single crystals. Using small angle neutron scattering, we have observed a square 
vortex lattice with the nearest-neighbors oriented at 45° from the Cu-O bond direction, which is 
consistent with theories based on the d-wave superconducting gap. However, the square symmetry 
persists down to unusually low magnetic fields. Moreover, the diffracted intensity from the vortex 
lattice is found to decrease rapidly with increasing magnetic field. 



It is a matter for debate whether hole-doped and 
electron-doped high-T c cuprate superconductors (HTSC) 
can be described within a unified physical picture |lj, 
2]. Indeed, electron-doped HTSC have markedly dif- 
ferent properties from hole-doped HTSC. For example, 
electron-doped materials have comparatively low values 
of the superconducting transition temperature, T Cl and 
much lower values of the upper critical field, B C 2- Fur- 
thermore, their normal-state resistivity varies as T 2 as 
expected for a Fermi-liquid [3|, |4| , and the presence of a 
pseudogap is still under discussion |5j]. Electron-doped 
HTSC also appear much closer to long-range antiferro- 
magnetic (AF) order, which can in fact coexist with su- 
perconductivity [1 0. In hole-doped HTSC, the 
d-w&ve nature of the order parameter is well-established 
|lflj . However, the evidence for the symmetry of the su- 
perconducting gap in electron doped materials (which 
has important implications for the pairing mechanism Q) 
is somewhat contradictory. Earlier measurements of the 
penetration depth and tunneling experiments [l2T ] 
supported s-wave symmetry, whereas a d-wave super- 
conducting order parameter is indicated by more recent 
phase-sensitive E3 an d ARPES experiments [Til Il5| . 
The electron-doped superconductors are of particular in- 
terest in this respect, since they have a tetragonal struc- 
ture (rather than orthorhombic) and therefore should 
show pure d-wave behavior, unaffected by admixture of 
an s-wave component associated with the orthorhombic- 
ity H3. 

Recently, there has been considerable interest in the 
nature of the vortex lattice (VL) in unconventional su- 
perconductors. For instance, vortex cores in d-wave su- 
perconductors are pre dicted to have a distinctive four- 
fold structure [3 Hi El El El] This leads to the 
expectation that a square VL is formed at high mag- 
netic field, with the nearest-neighbor directions aligned 
with the nodes of the order parameter. These theo- 
retical predictions are consistent with the small-angle 



neutron scattering (SANS) observation of a transition 
from an Abrikosov-like hexagonal VL to a square VL in 
the hole-doped HTSC YBa 2 Cu 3 07 (YBCO) 21]. How- 
ever, similar measurements on overdoped La2- a: Sr :l: Cu04 
(LSCO) [2^ show a square lattice with nearest neighbors 
oriented at 45° to the nodes of the superconducting or- 
der parameter. The orientation of the nodal directions in 
the heavy-fermion superconductor CeCoIns is at present 
uncertain (2^, so it is unclear if the recent SANS obser- 
vations of the VL in this material |2J| confirm theoretical 
expectations. It is therefore of great interest to gain fur- 
ther information from another class of superconductors, 
such as the electron-doped HTSC. Moreover, these com- 
pounds do not have twin planes present in orthorhombic 
systems such as YBCO and LSCO. These planar defects 
with a suppressed order parameter are capable of pin- 
ning the VL [2j|. Finally, the low values of the upper 
critical field (B c2 ~ 10 T compared to B c2 ~ 100 T in 
hole-doped HTSC) allow the investigation of the whole 
magnetic phase diagram. 

Only recently could large enough electron-doped crys- 
tals be produced, and to our knowledge, a SANS in- 
vestigation of their VL has not yet been published. 
We report here the first direct observation of a VL 
in Ndi.85Cen.i 5 Cu04 (NCCO). Our SANS experiments 
(see ref. [26| for a description of the technique) were 
performed on the instrument D22 at the Institut Laue 
Langevin, France, using neutrons with a wavelength X n = 
6 A- 20 A. Crystals of NCCO were grown in a mirror fur- 
nace and annealed as described in ref. |27| to give an on- 
set T c 25 K (AT C «3K, 10%-90% criterion). Samples 
with two shapes were investigated in applied magnetic 
fields up to 0.4 T. In high fields, a cylinder of 5 mm 
diameter, consisting of two nearly aligned crystals, was 
mounted in a cryostat with the magnetic field direction 
bisecting the two c-directions and at 7° to each of them. 
Single crystal plates of ^1.5 mm thickness were used at 
low fields: this required a long neutron wavelength and 
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the smaller thickness reduced the effects of neutron ab- 
sorption in the sample. They were mounted so that the 
c-direction was within ~ 2° of the applied field, which 
was approximately parallel to the incident beam. There 
were no significant differences in the VL diffraction pat- 
terns obtained from different samples in measurements 
at the same field. 

In Fig. H we show the VL diffraction patterns obtained 
at various magnetic fields. A background measured above 
T c has been subtracted because the VL signal is ex- 
tremely weak (~0.013% of the total detector counts at 
50 mT). For magnetic fields larger than 50 mT, the VL 
clearly has square coordination. The nearest-neighbor 
directions are parallel to the {110} crystal directions, 
which correspond to the nodes of the superconducting 
order parameter, in agreement with predictions for d- 
wave HTSC [3 Uji El E HJ Further evidence for a 
square VL is given by the positions of the Bragg spots in 
reciprocal space. For first-order diffraction, the relation- 
ship between the magnitude q of the wavevector and the 
magnetic field B is given by flux quantization: 

q = 27r \/-l-' W 
V ct$ 

where a is equal to 1 or a/3/2 for square and hexagonal 
lattices respectively, and $o is the flux quantum, q can 
be obtained by fitting the tangential sum of the neutron 
signal with a Gaussian, as shown in Fig. [21 As expected, 
the position of the peak shifts to higher q with increasing 
magnetic field (confirming the VL origin of the neutron 
signal), and the extracted values of a are consistent with 
a square VL at all fields measured (see inset to Fig. |2] for 
the lower fields). We have confirmed by magnetization 
measurements that the value of the trapped flux at low 
fields is within ~ 1% of the applied field. Hence, flux 
expulsion has an insignificant effect on the value of the 
q- vector. 

As shown in Fig. 2] the VL orientation is not perfect at 
any field. However, the orientation becomes rather more 
disordered as the field is decreased below 50 mT and 
at 20 mT the intensity distribution becomes ring-like. A 
similar distribution of intensity at low magnetic fields has 
been observed in LSCO 22], and was attributed to the 
superposition of diffraction patterns from various domain 
orientations of hexagonal coordination, since the value of 
a at low fields was consistent with that of an hexagonal 
VL. Moreover, by rotating the c-axis 10° away from the 
field direction, the degeneracy of the VL system could 
be reduced, and the hexagonal coordination of the VL 
in LSCO was confirmed. In NCCO, on the contrary, the 
values of a at low magnetic fields are still consistent with 
that of a square VL, and measurements at 20 mT with 
the c-axis rotated away from the field direction caused 
little change in the pattern. The effects of rotating of 
the field away from the c-axis at 50 mT are shown in 
Fig. |21 At 20° little change is observed, whereas at 30° 
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FIG. 1: SANS diffraction patterns measured at T w 2.5 K, 
after field cooling from 30 K in B = 20, 50, 100 and 200 mT. 
A background taken at T — 30 K has been subtracted. The 
{100} directions, corresponding to the Cu-O bond directions 
in the CUO2 planes, are oriented vertical/horizontal in the 
pictures. 
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FIG. 2: Tangential sum of the neutron signal vs wavevector q 
for B = 50 and 100 mT at T » 2.5 K. Inset: field dependence 
of the coordination-dependent quantity a (see Eq0 at low 
fields; the horizontal lines indicate the expected values for 
hexagonal and square VL. 

the pattern becomes rather disordered and the scattered 
intensity lies on an ellipse as expected from the uniaxial 
anisotropy of the crystal. However the average value of 
a is unchanged. It is surprising that the vortex lattice 
maintains square coordination without being aligned to 
the crystal lattice. 

In addition to the <i-wave scenario for square VL co- 
ordination, one should also consider the effects of Fermi 
surface anisotropy. An appropriate theory for large-K 
materials well below B C 2 is London theory with non- 
local corrections [i^, which has been extensively used 
to account for VL phase transitions in the borocarbides 

mils. 

If both these effects are present |3lJ , we would 
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FIG. 3: SANS diffraction patterns at B = 50 mT taken as in 
Fig.^but with the c-direction rotated about the vertical axis 
by 20° and 30° to the field direction. The growth direction 
and long axis of the sample lay 15° counterclockwise to the 
vertical and may be associated with the stronger intensity 
in these quadrants. The ellipses drawn have axial ratio of 
cos(30°). 
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FIG. 4: VL intensity measured at T ps 2.5 K as a function of 
magnetic field. The dotted line is the expected field depen- 
dence taking into account core effects and field-dependence of 
the penetration depth (see text). 



expect the strong angular variation of the d-wave gap to 
dominate over the usually smaller variation of the Fermi 
velocity. ARPES experiments together with band struc- 
ture calculations [32j indicate that NCCO and YBCO 
both have nearly isotropic hole-like Fermi surfaces with 
a slight four-fold distortion oriented so as to favor the 
observed square VL orientation. In overdoped LSCO, on 
the other hand, the Fermi surface is electron- like |33|, |3J] 
with a square shape oriented at 45° to that of the other 
two compounds. Moreover LSCO exhibits a pronounced 
anisotropy in the Fermi velocity. Both the shape of the 
Fermi surface and the Fermi velocity anisotropy would 
indicate (via nonlocal effects) the VL orientation actu- 
ally observed in LSCO 22]. 

However, both d-wave and nonlocal effects should only 
be important at fields which are substantial fraction of 
B c2 , when the inter- vortex spacing is comparable to the 
coherence length. For instance, <i-wave calculations in- 
dicate that the square symmetry has lower free energy 
than the hexagonal one at applied fields greater than 
0.15 B c2 w 1.2 T 0, or than B c2 /k w 0.4 T H3 (tak- 



ing the values of B c2 ~ 8 T [35J , penetration depth Xl ~ 
1250 A [3^| and coherence length £ w 60 A [33, there- 
fore k = ~ 20). Although the characteristic fields 
in electron-doped HTSC are generally lower than in the 
hole-doped materials, these estimates are too large to ex- 
plain our results. Hence, our observation of a square VL 
in NCCO down to very low magnetic fields is rather sur- 
prising, unless another source of anisotropy in the Cu02 
planes is present. One candidate for this is the Cu antifer- 
romagnetic correlations, whose characteristic wavevector 
7] coincides in direction with the VL reciprocal lattice 
wavevector. 

Another matter of great interest is the variation of the 
diffracted intensity with field. The intensity Ihk of a sin- 
gle (h, k) reflection (int egrated over the rocking curve of 
the VL) is given by [37| : 
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where <fi is the incident neutron flux, \x is the neutron 
magnetic moment, V is the sample volume and q^k is the 
(h, k) reciprocal lattice vector. Fhk is the "form factor" 
of the (h, k) reflection. It is a Fourier component of the 
spatial variation of the magnetic field, and in the London 
limit is related to the penetration length Al by: 



Fhk — 
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(3) 



For B c i <C B <C B c2 , the second term in the denom- 
inator is dominant and this gives Ihk °c </j~ fe oc B~?. 
In our NCCO samples, the rocking curves of the VL 
diffracted intensity were found to be fiat, over a range 
of ±4°, within experimental error. The expected width 
of the rocking curves can be estimated using Eqns. (0) 
and J2J and the value of Xl quoted in /iSR measurements 
■ At 50 mT we obtain a width larger than 10°, con- 
sistent with our experimental observations. We therefore 
measured the intensity as a function of field at a fixed 
sample angle. Assuming that the rocking curve width re- 
mains constant with field, the measured intensity is pro- 
portional to the integrated intensity. In contrast with the 
London prediction, we observe a strong field-dependence 
of the scattered intensity (see Fig. Eland Fig. 0}, which 
becomes immeasurably small above B w 0.4 T, which 
is well below B c2 . If represented by a power law, this 
variation has an exponent of about -2. (If the rocking 
curve becomes narrower with increasing field, then the 
field dependence of the VL intensity becomes even more 
rapid and more difficult to interpret. A broadening of the 
rocking curve as a function of increasing field would be 
consistent with increasing disorder as discussed below). 

We consider two possible intrinsic reasons for the fast 
decrease of the VL intensity with field. One is the effect 
of the finite size of the vortex core, which will lead to devi- 
ations from the London predictions at large values of the 
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wavevector q. This effect can be modeled by a calcula- 
tion using GL theory pjjj , giving an algebraic expression 
for the correction to the form factor: exp(—y/2q£) which 
may also be written as exp{—2(TtB/B c2 )2). In addition, 
there is expected to be a specifically d-wave effect - an 
increase of the penetration depth with field js^. /iSR 
measurements in YBCO indicate that the variation of A 
with B is linear at low fields (B < 2 T) |4(i|: 
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KB) 
A(0) 



1 + /3 



B 

B~ ' 



(4) 



where (3 is a temperature-dependent coefficient that re- 
mains finite at T—0 K, and Bq a characteristic field on 
the order of the thermodynamic critical field B c . If we 
assume that a similar linear dependence of X(B) is valid 
for NCCO, and include the core-size effects as well, we 
can model the variation of the VL intensity with field: 



hk ~ B- 



k ( 1 + P §~) 4ex P(-^B/B c2 )^) (5) 



For NCCO we used /3 = 7 ■ 1(T 2 (as for YBCO) and 
B ~ B c2 /V2k = 0.28 T. As can be seen in Fig. 0] 
we cannot represent the strong field-dependence of our 
experimental data by these intrinsic effects. It may be 
that an enhancement of AF correlations by the applied 
field is reducing the superfluid density. However, we 
find that even doubling (3 and the exponential prefac- 
tor does not reproduce our results. The rapid decrease 
of the VL intensity is more likely to be due to a tran- 
sition to a more disordered vortex system, as has been 
predicted theoretically 0J an d associated with a sec- 
ond peak in magnetization measurements (Z^. SANS 
experiments in the isotropic (K,Ba)Bi03 system 0] re- 
vealed a similar rapid loss of diffracted intensity near the 
second peak. In hole-doped Bi2.i5Sri.95CaCu208+a;, a 
strong decrease of the VL intensity with increasing field 
was attributed to a crossover from a 3D to a 2D vor- 
tex system 44] . In NCCO, however, such a dimensional 
crossover is expected to occur at much higher magnetic 
field B 2D =<I> /(js) 2 - 13 T (> B c2 l), where s ~ 6 A is 
the distance between Cu0 2 planes and 7 ~ 21 3] is the 
anisotropy. In this respect, NCCO seems to be similar 
to underdoped LSCO [45j, in which /nSR measurements 
have given clear evidence of a field-induced crossover to 
a more disordered, but still three-dimensional VL. 

In conclusion, we have made the first direct observa- 
tion of the VL in the electron-doped NCCO, which is the 
first tetragonal HTSC to be investigated by SANS. Con- 
trary to theoretical expectations, the VL remains square 
down to very small fractions of B c2 . In addition, we have 
observed an unusually fast decrease of the VL intensity 
with increasing magnetic field, which is probably due to 
a crossover to a more disordered vortex state. 

This work was supported by the Swiss National Sci- 
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